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Abstract. Weighted discrete Hilbert transforms (o„)„ i-» ( J2 n a nV n /(^j - In)) ^ from to 
are considered, where V = (7„) and A = (Xj) are disjoint sequences of points in the complex 
plane and v = (v n ) and w = (wj) are positive weight sequences. It is shown that if such a Hilbert 
transform is unitary, then T U A is a subset of a circle or a straight line, and a description of all 
unitary discrete Hilbert transforms is then given. A characterization of the orthogonal bases of 
reproducing kernels introduced by L. de Branges and D. Clark is implicit in these results: If a 
Hilbert space of complex-valued functions defined on a subset of C satisfies a few basic axioms 
and has more than one orthogonal basis of reproducing kernels, then these bases are all of Clark's 
type. 

1. Introduction 

If we are given two finite or infinite sequences of distinct points T = (7™) and A = (Aj) in C 
and a sequence of positive numbers v = (v n ), we may define the discrete Hilbert transform by 



To make sense of this, we assume that T and A, viewed as subsets of C, are disjoint. We also 
assume that A is a subset of the set 



If we now associate another weight sequence w = (wj) with A, we may ask: When is the discrete 
Hilbert transform H V (T, A) given by (1) a unitary transformation from if, to ^? 

The present note answers this question and shows how the solution translates into a general 
statement about orthogonal bases of reproducing kernels. Making a few minimal assumptions on 



(1) 





because we wish to define the discrete Hilbert transform for sequences (a n ) in 
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the underlying Hilbert space, we arrive at the following conclusion: There are no other orthogo- 
nal bases of reproducing kernels than those introduced and studied by L. de Branges [4] and D. 
Clark [3]. 

2. Localization of the sequences r and A 
Our starting point is the following localization result. 

Theorem 1. If the discrete Hilbert transform 

H V (T,A): tl^tl 
is unitary, then T U A is a subset of a circle or a straight line in C. 

Proof. In what follows, we let denote the vectors in the standard orthonormal basis for t\. 
Thus is the sequence for which the n-th entry is Vn 1//2 and all the other entries are 0. 
We fix an index m and observe that since T is a subset of (A, it?)*, the function 

Wj 



G(z) = (z- 7m ) £ 



(Aj -7m)(Aj - z) 

is well defined for z in Y. In fact, since H V (T, A) is assumed to be a unitary transformation, the 
basis vectors map into an orthonormal system in and therefore G vanishes on T. Thus 
we may write 

G(z) = G{z) - G( 7 ») = (z- 7n ) J2 rr -f! J " 7 "i ? 

j (Aj - 7m) (Aj - 7n)(Aj - z) 

where on the right-hand side we have just subtracted the respective series that define G(z) and 
Gr(7 n ). It follows that 



Z ~ln j 7m) (Aj - In) {Xj ~ z) 

and this function vanishes for z in T \ { r f n }. Since H V (T, A) is assumed to be unitary, the vectors 
H V (T, A)e^ constitute an orthonormal basis for i 2 w , and therefore the sequence 

Aj - 7 m 1 
Aj - 7m" Aj - 7fl 

is a multiple of the sequence (l/(Aj — tw))j. Thus the complex numbers (Aj — 7m) 2 /(Aj — 7 n )' 
have the same argument for all j, and so 
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for j ^ I and m ^ n. In other words, the cross ratio of the four complex numbers \j, Xi, j n , r ) m 
is real. This can only happen if the points lie on the same circle or straight line. 

After having applied this argument to four arbitrary points, say Ai, A 2 , 71, 72, we see that in 
fact every point from T U A lie on the circle or straight line determined by the four initial points, 
because we may apply the same argument to any given point in T U A along with three of the 
points Ai, A 2 , 71, 72- □ 



For a given sequence T being a subset of a circle or straight line and an associated weight 
sequence v, we wish to describe those pairs A and w such that H V (T, A) : — > t 2 w is a unitary 
transformation. To begin with, we require the admissibility condition 



which is now a necessary and sufficient condition for (T, v)* to be nonempty; we will say that 
(r, v) is an admissible pair whenever (2) holds. 

We will assume that T is a subset of the real line. The case when T is a subset of a circle is 
completely analogous, as will be briefly commented on at the end of this section. We set 



and observe that ip is well-defined on (T,v)* because the series in (3) converges absolutely for 
z in (r, v)*. We also note that (p is a Herglotz function in the upper half-plane ([2], Chapter 9 ). 
In particular, it is a holomorphic function whose imaginary part is positive. A general Herglotz 
function ip in the upper half -plane can be written as 



where b is a real constant, c a positive constant, and [i a nonnegative measure on the real line 
such that 



We will say that ip is a purely atomic Herglotz function if c = and \x is a purely atomic measure; 
our function </? is thus an example of a purely atomic Herglotz function. 
For every real number a, we set 



3. The unitary transformations associated with r and v 



(2) 




(3) 






A(a) = {A G (r»* : <p(\) = a}. 
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We observe that 

v n (z — w) . . . . 

(4) J2i — v r ^"^. 

which implies that the sequences (1/(A — 7 n ))n with A in A(a) constitute an orthogonal set in 
£ 2 . This means that A(a) is at most a countable set, so that we may associate with A(a) a weight 
sequence w(a) = (wj), where 

(5) Wj =(j2 



for Aj in A(a). It is implicit in our arguments that if H V {Y, A) : £ 2 — y £ 2 W is a unitary transforma- 
tion, then A = A(a) and w = w(a) for some real number a. 
We will now prove the following theorem. 

Theorem 2. Let (Y, v) be an admissible pair with V a subset of the real line, and let a be a real 
number. Then the discrete Hilbert transform 

is unitary if and only if (a — ip(z))~ l is a purely atomic Herglotz function. 

Proof. In this proof, we will again use the standard orfhonormal basis vectors in l 2 ; we will 
denote the corresponding basis vectors in (. 2 W ^ by We will use the notation || • ||„ and || • \\ w 
for the respective norms in £ 2 and t 2 w . 

It is clear that the adjoint transformation to H V (T, A(a)) is again a discrete Hilbert transform. 
In fact, since T and A(a) are sequences of real numbers, we have H*(T, A(a)) = —H w t a \ (A, Y), 
where 

H w{a) (A(a),Y):£ 2 w{a) ^£ 2 v . 

Therefore, H V (Y, A(a)) is unitary if and only if both H V (Y, A(a)) and H w ( a \(A(a), Y) are iso- 
metric. Hence it suffices to check whether (H V (Y , A(a))e^) and (H w r a -\(A(a),T)f^) are or- 
thonormal sequences in respectively £ 2 W ^ and £ 2 . 

The orthogonality of the vectors H w ( a \(A(a), Y)f^ in £ 2 has already been verified (see (4)); 
it is just a consequence of the definition of A(a). Likewise, we have automatically 

ll^(a)(A(a),r)/W||2 = X;i^^ = l- 

So our task is to show that (H V (Y, A)e^) is an orthonormal sequence in £ 2 (A(a), w(a)) if and 
only if (a — (p(z))" 1 is a purely atomic Herglotz function. 



UNITARY DISCRETE HILBERT TRANSFORMS 



5 



We first assume that (a — (p(z)) 1 is indeed a purely atomic Herglotz function. It suffices to 
show that there is a real constant b such that 



1 , ^ / 1 A 



■W: 



a — ip(z) ' 3 \Xj — z 1 + XjJ ' 

where Xj are the points in A(a) and Wj are as in (5). Indeed, by symmetry, it will then follow 
that the numbers 7„ are solutions to the equation 



w 



1 A, 



J{ Xj-z 1 + X) 



so that the arguments already employed for the vectors H v) f a \(A(a) , F)fW apply similarly to the 
vectors H V (T, A(a))e (n) . 

We start from the representation (6), with no a priori assumption on the points Xj and the 
nonnegative numbers wj except the admissibility condition 



< oo; 



+ A 



our goal is to prove that the Xj are in A (a) and that the Wj are given by (5). We first observe that 
if we set z — Xj + iy, then we get, by restricting to imaginary parts, 



" " (Aj - In) 2 + 2/ 2 



whence 



In other words, the points Xj belong to (T,v)*. We now multiply each side of (6) by z — Xj and 
take the limit when z — Xj + iy and y — > + ; since Xj is in (r, v)* and <p(Xj) = a, this gives (5). 

Suppose, on the other hand, that (a — Lp(z))^ 1 is not a purely atomic Herglotz function and 
that the vectors H V (T, A(a))e^ constitute an orthonormal system in £^ a y We will show that 
this leads to a contradiction. To begin with, our assumption on (a — ip^z))" 1 implies that 

(7) = b + ^ w > &FI " TTAj) + C2 + £ (rb " TTt 5 ) *««• 

with /i a spectral measure such that /i({Aj}) = for every j and not both c = and ji = 0; the 
fact that the u>j are given by (5) can be proved as in the first part of the proof. 
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We now argue in the same way as above, reversing the roles of T and A (at). This means that 
we first show, by again restricting to imaginary parts, that 

v Wj + r d ^ < 

V (7n-A;) 2 + i-oo(7n-t) 2 " » 

for every n. We infer from this that both the sum and the integral on the right-hand side of (7) 
converge absolutely for z = 7 n . Indeed, the right-hand side of (7) vanishes for z = 7 n , and so if 
we put z = 7„ + i5 in (7), divide each side by iy, and let y tend to 0, we get 

w, f°° d/i(t) 



(in - Aj) 2 J.oo (7„ - t) 



Since we should have \\H V (T, A(a))e^ n ^ = 1, we have reached a contradiction unless \x = 0. 
On the other hand, if /i = and c > 0, then we also reach a contradiction because the condition 
for orthogonality of the vectors H V (T, A(a))e^ becomes 

W-^ ^-)=o 

■ V Ira \ 7n A j y 

for m 7^ n, and this is inconsistent with the right-hand side of (7) being whenever z = j n . □ 

A few remarks are in order. First, it should be noted that we may have (T, v)* fl R = even 
if (r, v) is an admissible pair. The following is an example. Pick a sequence of distinct prime 
numbers p t such that 

E-l/2 . 
Pi < °°- 

I 

SetT = \J lP ^Z\Z, and equip T with the weight sequence v obtained by placing a weight of 
magnitude p ; 2 at every point of the sequence p^l, \ Z. 

On the other hand, if T is a discrete subset of the real line, then H V (T, A(a)) : il —> £ 2 W ^ is 
unitary for every a with one possible exception: It fails to be unitary when 

^n7n 



f n < oo and a = — 



+ 7 2 ' 

1 in 

This statement follows almost immediately from Theorem 2. We get the exceptional case because 
the constant c in the representation (7) is obtained as 

1 



lim 



y^-oa iy(a — <p(iy)) 
If r is a subset of the unit circle, then the potential (3) should be replaced by 

(8) <p(z) = ~2^v n ——- 

& in 6 
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the analysis goes through in the same way, and we obtain a statement completely analogous to 
Theorem 2. Note, however, that for discrete sets T on the unit circle, there will be no excep- 
tional value for a because there is no linear term 'cz' in the general representation of a Herglotz 
function. Indeed, a Herglotz function if) in the unit disk is of the form 



where b is a real constant and [i a nonnegative measure on the circle. 

Finally, as will be seen in the last section of this paper, the unitary transformations obtained 
from Theorem 2 (and its counterpart for the unit circle) correspond precisely to Clark's orthonor- 
mal bases [3]. From this point of view, Theorem 2 is essentially a reformulation of Clark's 
theorem. 



Let be a Hilbert space of complex- valued functions defined on some set f2 in C. To begin 
with, we assume that ffl satisfies the following two axioms: 

(Al) If / is in ^ and /(A) = for some point A in Q, then we may write /( z) — f (z)(z— A) 

with /o a function also belonging to Jtf. 
(A2) J4? has a reproducing kernel k\ at every point A in 

We wish to describe those spaces Jf? which admit orthogonal bases of reproducing kernels. To 
avoid trivialities, we assume that the dimension of Jf? is at least 2. We note that this family 
of spaces is part of the much larger family of spaces Jif that admit Riesz bases of normalized 
reproducing kernels. Since each space of the latter kind can be equipped with an equivalent norm 
such that one of the Riesz bases becomes an orthonormal basis, the question of interest is when 
a space Jtf? has more than one orthogonal basis of reproducing kernels. It is therefore reasonable 
to introduce a third axiom: 

(A3) There exists a sequence of distinct points T = (j n ) in f2 such that the sequence of 
normalized reproducing kernels (Ac 7n /||/c 7 „ || #?>) constitutes a Riesz basis for M '. In 
addition, there is at least one point A in \ T for which n x ^ 0. 

The Riesz basis (Ac 7 „/||« 7n ||jr) has a biorthogonal basis, which we will call (g n ). Thus 
9n{lm) — whenever m ^ n. We fix an index n and set E(z) = g no (z)(z — 7 no ). It fol- 
lows from axiom (Al) that E{z)/(z 7n) = 9n i z ) {in 7no 

)9n {z)/{z - j n ) also belongs 

to J4f. We use the suggestive notation E'(j n ) for the value of this function at 7„. We have 
E'iln) 7^ because otherwise E(z) /(z — 7„) would be identically 0, which can only happen if 
all functions in Jf? vanish at every point in f2 \ T; this would contradict the last part of (A3). By 
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uniqueness of the biorthogonal sequence (g n ), we now have 

E(z) 

9n[Z) E'{ ln ){z- ln ) 
for every n. We call E, which is unique up to a multiplicative constant, the generating function 
for T = (jn). We may assume that E does not vanish at any point A in \ T, because then 
E(z) /{z — A) would be a vector in ffi vanishing at every point in T. Hence E(z) /(z — X) would 
be identically 0, which again would be in contradiction with the second part of (A3). 

The sequence g n is also a Riesz basis for Jjf, and therefore every vector h in Jt? can be written 

as 

E(z) 



(9) Kz) = ^2h(-y n ) — 



E'{ ln ){z- ln ) 
where the sum converges with respect to the norm of Jtf? and 



Win)? 



* ~ |U "m2 < °°' 

■7n\\je 



since point evaluation at every point z is a bounded linear functional, (9) also converges pointwise 
in f2. By the assumption that h i— > (^(7n)/||^7 n ||jr) is a bijective map from ^ to £ 2 , this means 
that 



|S'(7«)| 2 k-7n| 2 

whenever z is in f2 \ T. We see that the generating function E appears as a common factor in (9). 
Since E(z) ^ for z in f2 \ T, the function E can be divided out. 
We set 

Ik- II 2 



l^'(7n)| 2 

and observe that since f2 \ T is assumed to be nonempty, (10) implies that 

Vn 
+ |7n 



^-^ 1 + r/. 1 



We may now change our viewpoint: Given a sequence of distinct complex numbers T = (j n ) 
and a weight sequence v = (y n ) that satisfy the admissibility condition (11), we introduce the 
space Jf?(T, v) consisting of all functions 



m = £ 



for which 



2 _ \ ^ I \2 



UNITARY DISCRETE HILBERT TRANSFORMS 



9 



assuming that the set (T, v)* is nonempty. The reproducing kernel of j4f(T, w) at a point z in 

(I»* is 

MC)= ?(^-7^)(C-7n)- 

If A = (Xj) is a sequence in (r, v)*, then we associate with it the weight sequence w = (wj), 
where 



|Aj - 7 n | 2 

Consequently, (kx j /\\k\ j ||jf(ryy)) is an orthonormal basis for J4?(T, v) if and only if H V (T, A) : 
£l — >■ t 2 w is a unitary transformation. Thus from the two previous sections we conclude: 

If the space J^(T,v) has an orthogonal basis of reproducing kernels, then T is a subset of 
a straight line or a circle. Moreover, when T is a subset of the real line, the orthogonal bases 
of reproducing kernels for J^f(T, w) are obtained from the unitary transformations described by 
Theorem 2; an analogous result holds when L is a subset of the unit circle. 

5. Relation to Clark's bases 

We are now finally prepared to point out the correspondence between our description of unitary 
discrete Hilbert transforms and the bases studied by de Branges [4] and Clark [3]. We restrict to 
Clark's bases; the only difference between the two cases is that Clark considered the case of the 
unit circle while de Branges worked on the real line with, in our terminology, \j n \ — > oo. 

Suppose Lf is of the form (8) with T = (j n ) a sequence of distinct points on the unit circle. 
Then the function 

v ; <p(z)+i 

is an inner function in the open unit disk D. We associate with / the so-called model space 
K] = H 2 IH 2 , which is the orthogonal complement to the shift-invariant subspace IH 2 of the 
Hardy space H 2 of the unit disk. Since 1/(1 — (z) is the reproducing kernel for H 2 at a point ( 
in D, the reproducing kernel for K 2 at the same point Q is 



1 " HOW 



l-Cz 

This formula extends to each point on the unit circle at which every function in K] has a radial 
limit whose modulus is bounded by a constant times the H 2 norm of the function. 
A computation shows that 

■ ! + /(*) ( , 
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which according to Clark's theorem means that the reproducing kernels 

!-/(;) 

constitute an orthogonal basis for Kj. In fact, Clark's theorem says that if (3 is a point on the unit 
circle and the spectral measure of the Herglotz function 

.(3 + I(z) 

(Pf)(z) = %- 



P - I(z) 

is purely atomic, then the reproducing kernels associated with the spectrum of (fp also constitute 
an orthogonal basis for K]. The spectral measures generated in this way correspond precisely to 
the spectral measures of the functions \/{a — <p(z)) with a any real number. 

Having observed this correspondence, we conclude that a Hilbert space Ji? of the type con- 
sidered in the previous section can have more than one orthogonal basis of reproducing kernels 
only if J^f is, up to trivial modifications, a model space Kj either in the unit disk or in the upper 
half-plane. 1 

An additional wonder, which can be seen from Clark's theorem or indeed by a straightforward 
computation, is that the norm in Jf? can always be computed as an L 2 integral over a circle or a 
straight line. 
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It may be objected that we fall short of characterizing all Clark bases because there exist model subspaces 
that do not satisfy axiom (Al) above. However, axiom (Al) can be relaxed in an appropriate manner so that this 
shortcoming is removed; we refer to our recent paper [ I ] for details. 
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